Josephson Frequency Singularity in the Noise of 
Normal Metal — Superconductor Junctions 
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A singularity at tlie Joseplison frequency in the noise spectral density of a disordered normal 
metal - superconductor junction is predicted for bias voltages below the superconducting gap. The 
non- stationary Aharonov-Bohm effect, recently introduced for normal metals H, is proposed as a 
tool for detecting this singularity. In the presence of a harmonic external field, the derivative of the 
noise with respect to the voltage bias reveals jumps when the applied frequency is commensurate 
with the Josephson frequency associated with this bias. The height of these jumps is non-monotonic 
in the amplitude of the periodic field. The superconducting flux quantum enters this dependence. 
Additional singularities in the frequency dependent noise are predicted above gap. 
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The observation of the Josephson effect Q typically re- 
quires two superconductors in contact. Here, we present 
two situations where a Josephson frequency can be ob- 
served in a normal metal - superconductor (NS) junc- 
tion. Although the Josephson frequency does not mani- 
fest itself in the average current through the NS junction, 
the noise characteristics bears clearly its signature. Su- 
perconducting features have been predicted in NS junc- 
tions: the doubling of the shot noise Q and the crossover 
from thermal noise to excess noise at iksT — {2e)V {V 
applied bias) are examples P,0|. Recent experiments 
in SNS junctions seem to be in qualitative agree- 
ment with this crossover. First we present a general 
framework to investigate the finite frequency noise of 
NS junctions, second, we propose an equivalent zero fre- 
quency, constant bias measurement in the presence of 
an harmonic perturbation to detect this analog of the 
Josephson effect. In the first case, the singularity at the 
Josephson frequency originates from the time oscillations 
of the current-current correlation function. In the sec- 
ond scheme, the superposition of an alternating field to 
the bias voltage leads to steps/singularities in the noise 
derivatives as a function of DC voltage. The strength 
of these singularities is non-monotonic with the ampli- 
tude of the harmonic perturbation. This so called "non- 
stationary Aharonov-Bohm effect" has been predicted 
PI and observed experimentally 0] in normal mesoscopic 
samples. Similarly, the singular behavior of finite fre- 
quency noise is understood both theoretically and exper- 
imentally |n,p[. For an NS junction with a bias smaller 
than the gap, the physical quantities can in principle be 
obtained from the normal metal results by replacing ev- 
erywhere the electron charge by the charge of a Cooper 
pair. In both the conductance and the excess noise of NS 
junctions, this doubling of the electron charge has to be 
divided by two in the prefactors, in order to account for 
the lack of spin degeneracy. However, no detailed calcu- 
lation of these effects for finite frequencies is available. 



and the generalization to above gap voltages cannot be 
obtained from an effective carrier charge. 

A general framework is provided for both the station- 
ary and non-stationary problems. Consider the coherent 
normal metal - superconductor junction in Fig. n^. A 
steplike dependence of the gap is assumed at the bound- 
ary, and the one channel states for electrons and holes 
incident from the normal side are specified by a scat- 
tering matrix with elements Sap {a,(3 = e,h). A time 
dependent phase <i>(t) is accumulated on the normal side 
by electrons impinging on the superconductor (incident 
holes will bear the opposite phase). Normally reflected 
particles will not accumulate this external phase. No in- 
elastic effects other than this perturbation are assumed 
close to the NS boundary. The probability for Andreev 
reflection [B is Ra = |s;ieP- In general, a time depen- 
dent formulation ||l3| of the Bogolubov-de Gennes (BdG) 
equations is needed. Here, we assume that the time de- 
pendence is adiabatic, so that appropriate elements of the 
scattering matrix are simply multiplied by exp[±2i$(t)] 
in the BdG equations ||ll|. The electron and hole wave 
functions are given by: 
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where v± — hk±/m is the velocity of waves with hk± — 
•\/2m(/^5 ± e). In what follows, we neglect the difference 
between the hole and electron wave numbers in Eq. (|l|) 
and (g), assuming the chemical potential fis to be large. 
Performing the Bogolubov transformation on the cur- 
rent operator, the statistical average of the current- 
current correlator is computed: 
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(updu*^)t^{u*pdUa)t2 + {Vf3dvl)t^{vpdVa)t2 
+ {Ui3dU*^)ti{v*^dVa)t2 + ivi3dVa)ti{u}dUa)t 

(w^aOti {UadVp)t2 + iupdVa)t2 

{UadVp)ti+ {updVa)ti {u*pdv*a)t2\ ■, (3) 



where udv — udxV — vdxU, and (( )) means that the 
square of the average current has been subtracted. The 
time dependence indicated with the indices ii and ^2 is 
exphcited in the electron and hole wave functions of Eq. 
(§ and (|). /e(£) = /(e - eV) {h{e) = 1 - f{-e - eV)) 
denote the electron (hole) distribution function and / is 
the Fermi-Dirac function. While only the terms con- 
taining the product /q(1 — fp) contribute to the low 
frequency noise |12| , |l3| , terms proportional to fafp and 
(1 — /a)(l — ffj), which involve matrix elements between 
states which differ by two quasiparticles, contribute to 
the finite frequency noise and to the non-stationary cal- 
culation below. The current-current correlator of Eq. 
(^ would have a spatial dependence if the exact wave 
vector of electrons and holes was included in this finite 
frequency calculation jlj]. However, this spatial depen- 
dence is only relevant at frequencies large compared to 
the inverse time of flight through the sample, which are 
not considered here 
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FIG. 1. a) Zero order spectral density of noise for the 
non-stationary Josephson effect (top) with a finite linewidth, 
and for the NS junction (bottom), b) The NS boundary. The 
light shaded region determines where the electron/hole wave 
function may accumulate phase: schematic description of the 
two scattering processes, Andreev and normal reflection. 

In order to describe non-stationary situations, we in- 
troduce the double Fourier transform of the current- 
current correlator: 



sini,n2)= / /rfiidt2e^("i*i+"^*^)((/(ti)/(i2))). (4) 



If the translational invariance in time is broken either 
by an external alternating field (with frequency 51) or 
spontaneously, like in the non-stationary Josephson effect 
(NJE) [y, this double Fourier transform can be written 
as: 

3(^1,^2)^ Y^ 2TTS{ni + n2-mn)s^"'\n2) . (5) 
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In both cases the zero harmonic, which is propor- 
tional to S{Qi + Q2) is the most standard quantity to 
study. In the presence of this time invariance this is 
the only term which is present. For the above men- 
tioned NJE, the zeroth order spectral density S^^'^luj), 
defined in Eq. (^ is proportional to the spectral func- 
tion line-shape D{lu), shifted by the Josephson frequency: 
5(°)(w) = D{uj - 2eV/h) which is depicted in Fig. 1 a). 

We first consider the finite frequency noise in the pres- 
ence of a constant bias. Because of translational invari- 
ance in time, the Fourier transform reduces to: 



^(r^i, r^a) = 27r^(f7i + n2)S{n2) . 
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where S{il,2) is the "usual" frequency dependent noise. 
In what follows, we assume that the energy dependence 
of the scattering matrix coefficients can be neglected for 
biases small compared to the gap. At arbitrary temper- 
ature and frequency below the superconducting gap, the 
noise is: 
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+f{e + eV + nuj){l-f{e~eV)) . (8) 



Note that the (constant) phase factor $ does not appear 
in the above because the geometry is open. The first 
term in Eq. (0) is important when thermal fluctuations 
are present. At zero temperature and for positive fre- 
quencies, only the second term contributes, this gives: 
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with the convention eV > 0. The striking feature in 
Eq. (pp is the singularity at the Josephson frequency 
2eV/h. This frequency scale appears in this normal su- 
perconducting geometry because below the gap, the only 
available charge transfer process involves the conversion 



of electrons into holes. To gain a better understanding 
of this singularity, remember that for a junction between 
two superconductors, the order parameter on each side 
oscillates as exp[—i2/X5j 2^/71] with ^Si and /is2 the chem- 
ical potentials on each side, leading to a current oscil- 
lation with frequency 2{fj,s2 ~ A'Si)/^- For a junction 
between two normal metals, where the electrons wave 
functions on each side oscillate like exp[— z/xi^2i/?i], a sin- 
gularity in noise at {fi2 — Mi)/^ has been pointed out 
m. In the case of an NS junction, the electron and hole 
wave functions on the normal side have a time depen- 
dence exp[— i(/i5 ± eV)t/h] (where fis is the chemical 
potential of the superconductor), leading to the singular 
behavior in Eq. (|9|) . A similar but weaker singularity has 
been pointed out in the fractional quantum Hall regime 
(FQHE) im, at the "Josephson" frequency e*V/h with 
e*/e the electron filling factor. 

At cj = 0, we recover the doubled shot noise of NS 
junctions although the current-current correlation func- 
tion in Eq. (jj) has not been symmetrized here. At zero 
frequency, the choice of considering {{I{ti)I{t2))) or its 
symmetrized analog does not matter. At finite frequency, 
depending on the measurement procedure, both correla- 
tors may occur [|6[. Here, the time dependence of the 
symmetrized correlator at i > h/eV is: 

((/W/(0)+/(0)/(t)))^^i..(l-i..).^-'(^^*/^) 



t2 



(10) 



Up to our knowledge experimental techniques do not al- 
low for the direct observations of these oscillations in 
time. Nevertheless, finite frequency measurements are 
possible B. 

We now turn to a situation where a time dependent 
vector potential is applied near the boundary (Fig. 1). 
Because low frequency measurements are more accessible 
than finite frequency ones, we suggest that a zero fre- 
quency noise analysis of a sinusoidally perturbed system 
is a more straightforward tool to investigate the pres- 
ence of the Josephson frequency in this NS system. In 
addition, this setup gives us an extra opportunity to ob- 
serve a superconducting behavior, i.e. a doubling of the 
electron charge, as shown below. The phase is chosen 
to be a periodic function of time $(i) = ^a sin(rii) with 
<i>a = 27rJ ^ dxAx/4>o where (p^ — hc/e is the normal 
flux quantum, [cci, X2] defines the interval where the vec- 
tor potential is confined (Fig. M. The effect of this 
perturbation on the average current is straightforward 
in the limit where Ra depends weakly on the energy: 
it brings a periodic modulation of the current A/ = 
(4e2//i)i?^[?ir2/e]$(jCos(51t). In the current-current cor- 
relations, this modulation leads to a non-monotonic effect 
as a function of phase, in contrast with the electromotive 
force action on the current. Note that no closed topology 
is imposed, in contrast to the usual AB effect. Neverthe- 



less, the phase is periodically modulated, and yields a 
non-zero contribution in the noise harmonics despite the 
open geometry. 

The low frequency noise is computed by first perform- 
ing the time integrals in Eq. (0) from the explicit time 
dependence of the current matrix elements. Using the 
generating function of the Bcsscl functions Jn, one ob- 
tains: 
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The phase 2$ (Fig. nl) accumulated in the Andreev pro- 
cess leads to a factor 2 in the argument of J^, which 
is reminiscent of the Cooper pair charge. The tempera- 
ture dependence in Eq. ( pTJ ) is specified by Fv{mh^) = 
{2eV — mhi^) coth[(2eF — mMl)/2kBT] determines how 
the steps in the noise derivative 
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are smeared with temperature. In Eq. (03) the sum 
over harmonics has a cutoff at M = YleV/Ml\ . In ex- 
periments B it is more convenient to characterize the 
non-monotonic dependence on voltage by taking the sec- 
ond derivative of the AB contribution to the noise. This 
is illustrated for two distinct temperatures in Fig. 2. 
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FIG. 2. Non-stationary AB effect: plot of d^S/dV^, ex- 
pressed in units of (8e*/7r7i^f2)_RA(l — Ra), as a function of 
2eV/hn, with the choice $„ = 3. For 2fcsr = 0.2hD. (full 
line) and for 2kBT — TiQ, (dashed line) 

For small temperatures Ti^ > 2kBT, one observes os- 
cillations as a function of 2eV/h.il, with a clustering of 
large amplitude peaks. For temperatures hil < 2kBT, al- 
though individual peaks can no longer be identified, clus- 
ters of of "large" steps in the noise derivative continue to 



give an average contribution to the non-stationary AB 
effect. This robustness enhances the likehhood of exper- 
imental observation, which is addressed below. 

In Ref. B, a small alternative microwave voltage 
Vi{t) — Vism{il,t) was superposed to the DC bias in- 
stead of applying locally a magnetic flux. If the drop 
in voltage occurs at the NS boundary then gauge invari- 
ance states that this is equivalent to the present proposal, 
with the substitution $a — eVi/hi}. If this drop is more 
extended, the analogy is restricted as the time of flight 
of electrons and holes through the relevant region is in- 
creased. Moreover, the experiments of Ref. ||] were per- 
formed in the diffusive regime, with a diffusion time td 
of the order of the fi"^. Although this seems to be at 
the border of applicability of the present approach, the 
assumptions made in the superconducting case are simi- 
lar, so one expects these structures to be rather robust. 
Note that for these predictions to be valid, the harmonic 
perturbation should not affect significantly the distribu- 
tion function of electrons in the reservoirs. In addition, 
for A ^ eV, a multichannel generalization of these re- 
sults is obtained straightforwardly with the substitution 
Ra{^ - Ra) = Z]ti^a„(1 — -Ra„), where Ra^ are the 
eigenvalues of the Andreev reflection matrix H,|l^,|l^,|l3[ . 

The present calculations are limited to the subgap 
regime, but can be readily extended to eV > A provided 
that one takes into account quasiparticle scattering states 
in the superconductor. In turn, a specific model for the 
NS boundary has to be chosen in order to specify the 
energy dependence of the scattering matrix coefficients. 
Calculations using the BTK pSJ model of a step pair po- 
tential with a delta function impurity potential will be 
discussed elsewhere in detail [|l^. Results indicate addi- 
tional singularities to the one at w = 2eV/h due to the 
presence of new scattering channels: electron-like quasi- 
particle transmission and hole-like (Andreev) quasipar- 
ticle transmission in the superconductor, which give rise 
to cusps at cj = {eV =F ^)/^- These occur because of 
the presence of a sharp edge in the density of states of 
the superconductor. Note that contrarily to the subgap 
regime, these striking features cannot be interpreted by 
the substitution of an effective charge. The above gap 
regime provides an additional justification for the sys- 
tematic study presented here: these frequency scales are 
tied to the energy intervals which are relevant for a given 
bias voltage, rather than an effective charge. Further in- 
creasing the bias A <^ eV ^ the superconducting singular- 
ity oj — 2eV/h weakens, and the dominant singularities 
{eV ± A)/?i « eV/h give the analog of the singularity of 
normal metals. 

The non-stationary Aharonov-Bohm effect provides an 
extra tool for the study of dynamical effects in low fre- 
quency noise. Because the effective charge 2e appears 
both the in the position of the steps in the noise deriva- 



tive and in the magnitude of the steps, this method of 
analysis could be applied to other quantum fluids with 
an effective charge which is different than e because of 
collective excitations. Shot noise measurements recently 
reported a direct measurement of the fractional charge 
pO| . The non-stationary Aharonov-Bohm effect could 
therefore also be applied for the investigation of frac- 
tional charge Josephson frequency signatures in the Hall 
effect. 

Part of G.B.L.'s work was done in Zurich, supported 
by the Swiss NSF. 
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